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ABSTRACT 

We present observations with the INTEGRAL/IBIS telescope of the wind nebula powered by the 
young pulsar B1509-58 and we discuss the spatial and spectral properties of the unpulsed emission in 
the 20 — 200 keV energy band. The source extension and orientation along the northwest-southeast 
axis corresponds to the jet emission seen at keV and TeV energies. The hard X-ray spectrum is 
consistent with the earlier Beppo-SAX measurements. It follows a power law with a photon index 
a = —2.12 ± 0.05 up to 160 keV. A possible break at this energy is found at the 2.9<r confidence 
level. The 0.1-100 keV data are consistent with synchrotron aging of pairs in the jet and yield a 
magnetic field strength of 22-33 /zG for a bulk velocity of 0.3-0. 5c. The synchrotron cut-off energy 
thus corresponds to a maximum electron energy of 400-730 TeV. 

Subject headings: acceleration of particles - gamma rays: observations - shock waves - pulsars: indi- 
vidual(PSR B1509) 



1. INTRODUCTION 

The young radio pulsar B 1509-58 is associated 
with the supernova remnant MSH 15-5,2 (G320 4-1.2) 
ijSeward fc Harnden 1981 iManchester et al. 1 9821. Its 
spin parameters (a period of 150 ms and period deriva- 
tive of 1.5 x 10~ 12 s s~ x ) make it one of the youngest 
and most energetic pulsars known, with a characteristic 
age r ~ 1700 yr and a spin-down power of 1.8 x 10 37 
ergs s _1 . It presents o ne of the largest m agnetic field 
strength (1.5 x 10 13 G) l|Kaspi et al. 19941) recorded for 
an isolated pulsar. Recent X-ray and 7-ray images have 
revealed a very complex pulsar environment where one 
can study several manifestations of pulsar wind nebu- 
lae: the shocked wind in an equatorial flow, a powerful 
jet, and the interaction and confinement of these wind 
features within the surrounding remnant. 

ROSAT first detected the wind nebula in X 
rays, showing its 17'-long tail extending to the 
southeast along a faint rid ge of radio emissio n 
(|Greiveldinger fc Qgelman 19931 iTrussoni et al. 1 9961. 
later confirmed by BeppoSAX. The tail length decreases 
from 17' to 8.5' (at 10 % level of the central nebula inten - 
sity) between 1.6 keV and 10.0 keV l|Mineo et, al. 2001|) . 
High-resolution Chandra im ages revealed further details 
of the non-thermal structure l)Gaensler et al. 2002(1 : first, 
the axis of the elongated nebula (with a position an- 
gle of 150° ± 5° north through east) which may corre- 
spond to the pulsar spin axis; then two northern toroidal 
arcs of emission, 17" and 30" away from the pulsar, 
that may correspond to wisps in an ion-loaded equa- 
torial flow as in the Crab nebula; and a 4'-long bright 
collimated feature along the axis which is interpreted as 
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a jet with a velocity Vj P .t > 0.2c llGaensler et al. 2002t 
iTamura et al. 19961 iBrazier fc Becker 19971) ! The ap- 
parent absence of a counterjet sets a lower limit of 5 to 
the brightness ratio between the two sides (if two jets are 
powered) . Doppler boosting can account for this contrast 
provided that the jet bulk velocity and its inclination to 
the line of sight verify (3j et cos9j et ~ 0.28. The wind car- 
ries away at least 0.05% of the pulsar spin down power. 

The emission found by H.E.S.S between 280 GeV 
and 40 TeV presen ts a comparable morphology 
ijAharonian et al. 2 0051. The elongated nebula extends 
on both sides of the pulsar along the X-ray axis, over 
6' (1 a level) to the northwest and southeast direction. 
The 7-ray spectrum is well fitted by an ^- 2 - 27 ±0- 2 3±o.20 
power law and can be explained by TeV electrons up- 
scattering the cosmological microwave background and 
the ambient interstellar radiation field. Fitting both the 
TeV emission and the BeppoSAX one gives a mean mag- 
netic field of 17 /zG which is in agreement with lower 
limits found for the j et in the equipartition assumption 
I DcLan£v_£kiiL_2QQ£ )■ At a distance of 5.2 ± 1.4 kpc 
i Gaensler et al. 1999(1 . the bright part of the X-ray jet, 
the longer southeast X-ray tail, and the coincident TeV 
one extend from 5 to 30 pc on the sky. 

We present observations of the PSR B1509-58 region 
obtained with the INTEGRAL/IBIS telescope operating 
in the hard X-ray range, describing the morphology and 
spectral properties of the emission and comparing them 
with the X-ray and TeV data. 

2. OBSERVATION AND DATA ANALYSIS 

PSR B1509-58 was observed with the IBIS coded aper- 
ture telescope onboard the INTEGRAL spacecraft. It 
consists of a dual detection layer operating between 
~ 15 keV and 10 MeV, and of a tungstene coded mask 
located 3.2m above the detector. The first and pix- 
elized detector plane, ISGRI, is composed of 128 by 
128 Cadmium-Telluride (CdTe) semiconductor detectors 
covering the energy range from ~ 15 keV to 1 MeV 
l|Lebrun et al. 20031) . Sky images are obtained by de- 
convolution of the ISGRI image by the mask pattern. 
The resulting angular resolution is 5.7' {la). 
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The exposure time of about 1.3 Ms corresponds to ded- 
icated observations centered on PSR B1509-58 in 2005, 
and of Galactic Plane Scan (GPS) observations from 
2003 to 2005. It consists of more than 700 science win- 
dows of about 30 mn each. The analysis was done with 
the standard INTEGRAL OSA software (version 5.1). 

2.1. Phase resolved analysis 

In order to study the pulsar wind from PSR B1509- 
58 and to reduce contamination by pulsed emission from 
within the light cylinder, we have analyzed the data from 
the off-pulse phase interval. To extract the light curve, 
the photon UTC arrival times were converted into the 
solar system barycentre using the pulsar position. Pho- 
tons were then folded into 20 phase bins using the ATNF 
ephemerides 7 . In each phase bin sky images were ob- 
tained from standard deconvolution of ISGRI detector 
maps. This method provides an automatically back- 
ground subtraction and the unpulsed emission is thus 
clearly visible. 
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Fig. 1. — The B1509-58 pulsar light-curve between 20 and 50 
keV. 



The resulting light curve is presented in Figure ^ A 
significant unpulsed component is clearly seen in the 
0.77 — 1.07 phase interval similar t o the ones found by 
BeppoSAX l)Gusumano et al. 20011) that will be used for 
the following imaging and spectral analyzes. A detailed 
study of the pulsed emission is deferred to another paper 
(Hermsen et al., in preparation). 

2.2. Imaging extended sources with the IBIS/ISGRI 
telescope 

As shown in Figure [21 the source detected with 
IBIS/ISGRI in the off-pulse emission (17-40 keV) ap- 
pears to be slightly extended. In the standard OSA 
software, the flux of a point-source is given by the 
peak height of the associated Point Spread Function 

7 |http : / /www, atnf.csiro.au/ research/pulsar/archive| 



ijGros et al. 2003j) . This is, however, incorrect for an 
extended source because of the flux dilution over a re- 
gion larger than the PSF width. A method has been 
developed for extracting the flux and its associated error 
of an extended source seen by IBIS and more generally 
by any coded mask telescope using the MURA pattern 
ijRenaud et al. 2006j) . This method converts standard 
flux images into intensity images (ie flux per sky pixel or 
per steradian). The total flux is estimated by summing 
intensities over the whole source extent. The dilution 
factor will be defined as the ratio between the real flux 
and the one given by the standard INTEGRAL software. 
This factor depends on the source size. 

3. RESULTS 

3.1. Imaging : an extended source around PSR 
B 1509- 58 

Maps of the off-pulse emission have been constructed 
in two energy ranges, at 17 — 40 and 40 — 100 keV. The 
significance map in the low-energy band is presented in 
Figure|21 A bright source is clearly visible at the position 
of PSR B1509-58. It is centered on RA= 228.50° ± 0.03° 
and DEC= -59.43° ±0.03° and it is fully consistent with 
the pulsar location. Figure shows the source radial 
profile along the NW-SE axis in the significance map, 
together with the profile expected from a point source. 
The source appears to be slightly extended and its ori- 
entation coincides with the morphology known at other 
wavelengths. 
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Fig. 2. — Significance image of the off-pulse emission in the 17-40 
keV band around PSR B1509-58. The coordinates are in J2000. 
The black contours out line the TeV emission as seen by H.E.S.S 
lAharonian et al. 20051) . 

With coded mask imaging, sky image pixels are not 
statistically independent, but two pixels separated by 
more than the PSF width (apsF = 5.7' ) can be consid- 
ered as independent. After subtracting the PSF in the 
significance map (see Figure 0} we obtain a confidence 
level of 4cr and 2.3er in the low and high-energy bands 



Integral observation of the PSR B1509-58 PWN 
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Fig. 3. — Smoothed profile of the significance map along the NW- 
SE axis between 17 keV and 40 keV. The red curve corresponds to 
the expected profile in a point source case. 



from the two residuals measured along the NW-SE axis. 
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Fig. 4. — Smoothed significance map excess after point source 
contribution removal between 17 keV and 40 keV. The white con- 
tours outline the TeV emission as seen by H.E.S.S. and the blue 
ones correspond to Chandra observation. Two excesses visible in 
the jet direction (NW-SE) are in global agreement with the H.E.S.S 
contours baring in mind the smearing due to the PSF. 

The source orientation and extension can be evaluated 
by fitting the intensity images with a 2D Gaussian with 
separate widths along each axis and a free rotation angle. 
To compare with a true point-source the same fit has 
been applied to the XTE J1550 source in the same field 
of view. 



In the low-energy band, XTE J1550 is well fitted by a 
symmetrical Gaussian of width <txte — 5.85' ± 0.07'. It 
is fully consistent with the instrument PSF apsF = 5.7'. 
For PSR B1509-58, the best fit yields an asymmet- 
ric Gaussian with a rotation angle of 155° ± 4° from 
north through east. Perpendicular to this axis, the 
source width of 5.90' ± 0.08' is equivalent to that of a 
point source. Along the NW-SE direction, the width 
o~nw~se = 7.98' ± 0.07' suggests that the source ex- 
tends over several arcminutes. The apparent source 
width results from the convolution of the PSF and 
the true source size. The latter is found to be a s — 

( a NW-SE ~ u psf) = 5-53' ± 0.07' along the major 
axis. In the energy band from 40 to 100 keV, the source 
is smaller with a s = 3.52'±0.07' along the main axis. All 
the quoted errors include the statistical uncertainty and 
a crude estimate of the impact of the map binning and 
correlation between adjacent pixels obtained by moving 
the gaussian location in the fit. 

The fact that the source extension and orientation 
compare with the morphology seen in X rays and 7 
rays, as well as the increased brightness to the southeast, 
strongly suggest that the unpulsed emission is dominated 
by the wind nebula. 

3.2. Spectroscopy : a possible spectral break near 160 

keV 

The dilution correction factors due to the source extent 
with respect to a point source flux are respectively 1.36 
and 1.19 in the low and high energy bands. Values at 
intermediate energies have been interpolated. The spec- 
trum of the off-pulse emission has been measured from 
17 keV to 200 keV using the standard spectral analysis, 
then divided by the selected phase width, and corrected 
for dilution as described above. It is presented in Fig- 
urcElt ogether with the former BeppoSAX measurements 
( Minc o et al. 2 0011. The INTEGRAL spectrum is well 
fitted by an £;- 2 - 12 ±0- 05 power law in agreement with the 
BeppoSAX and RXTE results (Q). The total flux from 
17 to 200 keV amounts to 0.30 ± 0.02 keV cm" 2 s" 1 . 
The upper limit found in the 130-200 keV bin suggests 
a spectral break near 160 keV. This upper limit indeed 
deviates at the 2.9 a level from the extrapolation of the 
power-law spectrum in this energy interval. The break 
is not visible in the pulsed emission and so cannot be 
attributed to the pulsar. It is consistent with the Beppo- 
SAX measurements at high energy. 

3.3. Discussion and conclusions 

The extent and orientation of the unpulsed emission 
seen by INTEGRAL/IBIS nicely match the morphology 
of the wind nebula and jet seen at keV energies with 
ROSAT, Chandra, and BeppoSAX, and at TeV energies 
with H.E.S.S. The source width is not consistent with 
pure DC emission from the pulsar itself, although a large 
contribution cannot be ruled out. There is a close agree- 
ment between BeppoSAX and INTEGRAL in spectra an 
source elongation. All circumstantial evidence suggests 
that the measured hard X-ray emission is dominated by 
jet emission. The elongated source appears to be slightly 
brighter to the south, with a flux ratio of ~ 1.4 between 
the SE and NW sides. Larger contrasts are found in X 
rays (~ 5) and in TeV 7 rays (~ 3.5), but the difference 
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■ BeppoSax spectrum 

■ Integral/IBIS spectrum 
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Fig. 5. — Energy flux spectrum of unpulsed emission from the 
ISGRI and BeppoSAX data. All error bars and upper limits are 
given at the 1 <r confidence level. 



B. The spectral cut-off possibly seen at 160 keV then 
yields maximum electron energies in the jet of 400 and 
730 TeV, respectively. 
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The field and maximum energy estimates given above 
would not apply if the high-energy part of the INTE- 
GRAL source is dominated by the toroidal part of the 
wind instead of the jet as assumed above, because the 
rapid increase of the field strength to equipartition val- 
ues in the post-shock flow would cool the particles more 
efficiently and the INTEGRAL data would not relate to 
the larger scale ROSAT and HESS observations. 

Following a leptonic scenario and fitting both syn- 
chrotron and inverse Compto n radiation, t he HESS data 
yield a mean field of \l\iG (jAharonian et al. 2005|) in 
reasonable agreement with the above range. Increasing 
B to the above values would result in 2 or 3 times fewer 
particles in the modelled flow, but the inverse Compton 
TeV brightness could still be accomodated owing to the 
large uncertainty in the actual infra-red radiation field in 
the jet vicinity. 



may be due to the poorer angular resolution of ISGRI 
which cannot separate the central wind components from 
the larger scale jet. In the following discussion we will 
explore the consequences of a common origin of the 0.1- 
160 keV and TeV emissions in the collimated wind flow 
or jet along the NW-SE axis. 

The relative energy fluxes recorded in X rays and TeV 
7 rays indicate that s ynchrotron losses dom inate over 
inverse Compton ones ijAharonian et al~2 005). Figure 
also suggests possible decrease of the jet app arent length, 
as seen by ROSAT (|Trussoni et al. 19961) . BeppoSAX 
ijMineo et al. 200"H) . and now IBIS. In a simple scenario 
with a linear jet with constant cross-section and velocity, 
thus with uniform magnetic field and density to conserve 
the magnetic and momentum fluxes, the maximum dis- 
tance Lj et reached by the flow before the pairs burn off 

their energy scales as Lj et oc E^ 1 ^ 2 ijPacholczvk 1970(1 : 



-3/2^ f^y 



-1/2 



IkeV 



pc 



L jet = (1.69xl0 4 ) jBt ain(0 jet )(^) 

(1) 

Figure H3 shows that the data are consistent with this 
simple dependence and that radiative aging is effec- 
tive above keV energies. This aging is also consistent 
with the spectral turn-over seen in Figure [5] The best 

Lj e t E-y fit to the data in Figure El gives a measure 
of the product f3j e tsin(9j e t 

■^-3/2 The - et ou j . flow has 

a typical velocity of 0.5c l|DeLanev et al. 2006(1 equiva- 
lent to the speed observed in several other wind nebu- 
lae (Crab Nebula (0.4c), Vela (0.3c-0.7c), Gil. 2-0.3 (0.8- 
1.4c)). It is consistent with MHD simulations of rela- 
tivistic jets (0.5c). A lower limit of 0.3c is inferred from 
the Doppler boosting constraint (3j e tCos( ~ 0.28. For 
Pj e t values of 0.3 and 0.5, we find field strengths of 22 
and 33 /iG, respectively. The large uncertainty in the 
pulsar distance implies an additional 27% uncertainty in 
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Fig. 6. — Estimates of the jet apparent length, Lj et sin(8j et ), as 
a function of energy. 

In summary, the comparison of the morphology of 
the wind nebula above 20 keV with lower-energy data, 
and the observation of a possible cut-off near 160 keV 
are consistent with a simple jet scenario with efficient 
synchrotron cooling in a mean magnetic field of order 
22 — 33/iG and outflow velocities of 0.3-0. 5c. The pres- 
ence of electrons with still extreme energies of 400-730 
TeV at several parsecs from the pulsar places a useful 
constraint on the unknown jet acceleration process. 



We gratefully thank the Australian Pulsar Timing 
Archive for making all the pulsar cpchcmcridcs available 
to us. 



Integral observation of the PSR B1509-58 PWN 



5 



REFERENCES 



Aharonian, F., et al. 2005, A&A, 435, L17 

Brazier, K. T. S., & Becker, W. 1997, MNRAS, 284, 335 

Cusumano, G., Mineo, T., Massaro, E., Nicastro, L., Trussoni, E., 

Massaglia, S., Hermsen, W., & Kuiper, L. 2001, A&A, 375, 397 
DeLaney, T., Gaensler, B. M., Arons, J., & Pivovaroff, M. J. 2006, 

ApJ, 640, 929 

Gaensler, B. M., Brazier, K. T. S., Manchester, R. N., Johnston, 

S., & Green, A. J. 1999, MNRAS, 305, 724 
Gaensler, B. M., Arons, J., Kaspi, V. M., Pivovaroff, M. J., Kawai, 

N., & Tamura, K. 2002, ApJ, 569, 878 
Greiveldinger, C, & Ogelman, H. B. 1993, Bulletin of the American 

Astronomical Society, 25, 1345 
Gros, A., ct al. 2003, A&A, 411, L179 

Kaspi, V. M., Manchester, R. N., Siegman, B., Johnston, S., & 
Lyne, A. G. 1994, ApJ, 422, L83 



Lebrun, F., et al. 2003, A&A, 411, L141 

Manchester, R. N., Tuohy, I. R., & Damico, N. 1982, ApJ, 262, L31 
Marsden, D., et al. 1997, ApJ, 491, L39 

Mineo, T., Cusumano, G., Maccarone, M. C, Massaglia, S., 

Massaro, E., & Trussoni, E. 2001, A&A, 380, 695 
Pacholczyk, A. G. 1970, Series of Books in Astronomy and 

Astrophysics, San Francisco: Freeman, 1970, 
Rcnaud, M., Lebrun, F., Terrier, R., et al. 2006, A&Ain press 
Seward, F. D., & Harnden, F. R. 1982, ApJ, 256, L45 
Tamura, K., Kawai, N., Yoshida, A., & Brinkmann, W. 1996, PASJ, 
48, L33 

Trussoni, E., Massaglia, S., Caucino, S., Brinkmann, W., & 
Aschenbach, B. 1996, A&A, 306, 581 



